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Abstract: We employ a basic formalism from convex analysis to show a simple 
relation between the entanglement of formation Ep and the conjugate func- 
tion E* of the entanglement function E(p) = S(Tiap)- We then consider the 
conjectured strong superadditivity of the entanglement of formation Ep{p) > 
Ep(pi) + Ep(pn), where pi and pn are the reductions of p to the different 
Hilbert space copies, and prove that it is equivalent with subadditivity of E* . 
Furthermore, we show that strong superadditivity would follow from multiplica- 
tivity of the maximal channel output purity for quantum filtering operations, 
when purity is measured by Schatten p-norms for p tending to 1. 



1. Introduction 

One of the central quantities in quantum information theory is the entanglement 
cost of a state, defined as the number of maximally entangled pairs (singlets) re- 
quired to prepare this state in an asymptotic way. Calculating the entanglement 
cost of a general mixed state as such is, with the present state of knowledge, a 
formidable task because one has to consider an infinite supply of singlets and 
construct a protocol using local or classical (LOCC) operations only, such that 
the resulting (infinite-dimensional) state approximates an infinite supply of the 
required state to arbitrary precision. Furthermore, the protocol must have max- 
imal yield, the number of states produced per singlet. The entanglement cost is 
the inverse of this yield. 

An important theoretical breakthrough was achieved in X, where the en- 
tanglement cost Ec was shown to be equal to the regularised entanglement of 
formation: Ec(p) = lim^oo Ep(p® n )/n. The entanglement of formation (EoF) 
(defined below in Q) is defined in a mathematical and non-operational way and 
is therefore much more amenable to calculation. Moreover, for 2-qubit mixed 
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states, a closed formula for the EoF exists [2]. Nevertheless, calculating the en- 
tanglement cost still requires calculations over infinite-dimensional states. For 
that reason one would hope for the additivity property to hold for the EoF: 
E F (pi ® p 2 ) =7 Ef{px) + E F (p 2 ), because then Ec = E F . Additivity of the 
EoF has been proven in specific instances 3 , 4 , 5|E1IZ|IH1 • Some of these additiv- 
ity results are sufficiently powerful to allow calculating the entanglement cost 
for certain classes of mixed states jSJEHZJ- The much sought-after general proof, 
however, remains elusive for the time being and, in fact, general additivity is 
still a conjecture. 

ft is very easy to show that the EoF is subadditive: 

E F (pi <g> p 2 ) < E F ( Pl ) + E F {p 2 ). 

Additivity would then follow from superadditivity: 

E F {pi®p 2 ) >TE F ( Pl ) + E F (p 2 ). 

In 0] a stronger property, which would imply (super) additivity, has been con- 
jectured for the EoF, namely strong superadditivity: 

E F { P ) >? E F ( PI ) + E F (p u ), 

where p is a general state over a duplicated Hilbert space and pj and pu are its 
reductions to the different copies of that space. 

In this paper we show that strong superadditivity of EoF is equivalent to 
subadditivity of a much simpler quantity, the so-called conjugate of the entan- 
glement functional E(p) = S'(Tr^p). We then exploit this equivalence to show 
that strong superadditivity would follow as a consequence of multiplicativity of 
the maximal output purity, measured by a Schatten norm, for quantum filtering 
operations (this quantity will also be defined in due course). 

The main results are stated in Theorems I and 2. To arrive at these results, 
we have made use of a basic formalism from convex analysis |?HllOj and we hope 
that our results will stimulate usage of this elegant theory in other areas of 
quantum information. 

2. Notations 

Let us first introduce the basic notations. Let S(p) denote the von Neumann 
entropy S(p) — — Tr phi p. For state vectors we will typically use lowercase Greek 
letters, tp, <f>, and assign the corresponding uppercase letter to the projector of 
that vector; e.g. & — \\p){ip\. For mixed states we will use lowercase Greek letters 
p, a, t. The identity matrix will be denoted by I. 

We shall denote the set of bounded Hermitian operators over the Hilbert 
space Ti. by B S (H), the set of non-negative elements in B S (TL) by B + (7i), and 
the (convex) set of all states (trace 1 positive operators) over 7i by S(H). 

We will frequently slim down expressions like max p6< s{. . .} to max p {. . .}. 
When the domain of, say, a maximisation over states is missing it will be im- 
plicitly understood that the whole of state space S(H) is meant. The abovemen- 
tioned naming convention for states and vectors will be adhered to exactly for 
that reason. 
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Any state p can be realised by an ensemble of pure states. An ensemble 
is specified by a set of pairs {(pi,tpi)}iLi, consisting of N state vectors ipi and 
associated statistical weights pi (with Pi >0 and J^. pi — 1). Here, N is called the 
cardinality of the ensemble. The entanglement of formation (EoF) of a bipartite 
state p (i.e., a state over the bi-partite Hilbert space Ha is defined by 



3. Convex Closures 

Admittedly, the definition of the EoF just mentioned is not very handy to work 
with. Not in the least because for generic states p the cardinality N of the 
optimal realising ensemble must be larger than i? 15 /4, where R is the rank of 
p |12| . This is one of the reasons why no really efficient numerical algorithms 
have been found yet to calculate the EoF Furthermore, the mere fact that 
the minimisation involves ensembles at all makes a theoretical study of the EoF 
rather difficult. One of the first attempts at proving additivity of EoF relied on 
the investigation of these optimal ensembles [Hj • 

The results in the present work depend on the following simple observation. 
The import of the definition (JIJ of the EoF, as has been shown in |14ll4") . is 
that the EoF is the convex closure (or convex roof, as it is called in [14| ) of the 
pure state entanglement function E{&) = S(TrA restricted to the set of pure 
states. This means that the epigraph of the EoF (being the set of points (p, x) 
in S(TL) x M with x > Ep(p)) on the complete state space S(TC) is the convex 
closure of the epigraph of the function E' defined over S(H), where 



This follows immediately from Cor. 17.1.5 of [Hj and the definition (JJJ. Note 
now that E is concave over its domain. There is, therefore, no need to explicitly 
exclude mixed states 1 , so Ep is the convex closure of E as well. 

In the following paragraphs we will apply the standard convex analytical 
formalism for convex closures to general bounded functions / whose domain is 
the convex set of states S(7i). We will denote the convex closure of / by /. One 
definition of the convex closure of / is 



agreeing, indeed, with the definition of the EoF. A less cumbersome formulation 
of the convex closure is based on Cor. 12.1.1 of 0, which states that the convex 

Of course, E(p) has no real physical significance for mixed states. Moreover, we must 
be careful to distinguish between the two possible definitions E(p) = S l (Tr J 4 p) and E'(p) = 
S(Tvb p)- On pure states, these two definitions yield the same value, but for mixed states this 
is not so anymore. 
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closure of a function / is the pointwise supremum of the collection of all afhnc 
functions on S(H) majorised by /. So, for all states p: 

f(p)= sup {TrpX : (Wip E H : Tr&X < /(#))}. (3) 

XeB 3 (H) 

The mentioned afhne functions are here the functions Tr^X, where X ranges 
over B S {TL) 2 . This dual formulation is then further simplified by defining an 
intermediate function /*: 

f*(X)= max Tr[pX] - f(p), (4) 
pes(H) 

the so-called conjugate function of / p . If / is continuous, then the conjugate 
function is just the Legendre transform of /. The conjugate function is convex 
in X, because it is a pointwise maximum of functions that are affine in X . The 
importance of the conjugate function is that the conjugate of the conjugate of / 
is the convex closure of /, / = /**, and the conjugate of the convex closure of 
/ is the conjugate of /, /* = /* (UJ, the remark just before its Theorem 12.2). 
Thus 

f(p) = max . Tr[pX}-f*(X) (5) 

A to ( ri) 

f*(X) = max Tr[pX] - f(p). (6) 
P eS{H) 

In other words, the conjugate and convex closure determine each other com- 
pletely. 

Because /* and / are convex functions, the optimal X and p in © and ©, 
respectively, both form convex sets (possibly singleton sets). Furthermore, there 
is a correspondence between the optimal X in (JSJ and the optimal p in ©. 

Proposition 1. (a) If X' is an optimal X for t in then (i) r is an optimal 
p for X' in and (ii) all members of an optimal realising ensemble for t are 
optimal p for X' in Op. (b) If p' is an optimal p for Y in then Y is an 



optimal X for p' in 

Proof. Statement (a)(i) is proven by inserting in and exploiting the 
premise that X 1 is an optimal X. This gives /(r) = Tr tX 1 — max p (Tr pX'-f(p)). 
Putting p = t yields an upper bound on the right-hand side because r is not 
necessarily optimal in the maximisation. However, the value of the bound we 
obtain is /(t), which happens to be equal to the left-hand side. Thus this choice 
really is an optimal one, proving optimality of r for X' in ©. 
Statement (b) is proven similarly, by inserting (0 in 

Considering statement (a)(ii), let {(pi,Ti)} be an optimal ensemble for r 
(with Pt > 0). Thus /( T ) = EiPif ln)- By assumption, /(r) = Tr tX' - 
f*(X'). Inserting Q and expanding unity as Y^iPi yields J2iPif( T i) = TrrX' — 
J2i Pi max p (Tr pX' — /(/?)). If we now replace p by n in the ?-th summation term 
we get an upper bound on the right-hand side, with equality only if all the t, 

2 For our purposes the Corollaries from |5] have to be restated with R" replaced by S(7i). 
This causes no problems if one extends the domain of / to the affine space of all trace 1 
Hermitian operators and defines f(x) = +oo for negative x. 
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are optimal p for X' . The bound is easily seen to be X^Pj/( r i)j which is ac- 
tually equal to the left-hand side. We find again that the bound is sharp, and 
optimality of the r, follows. □ 



4. Additivity 

These basic results will now prove to be a powerful tool for studying the additiv- 
ity issue of the EoF. Let Hi and Hn be two copies of the Hilbert space Ha®Hb, 
and define H = Hi^Hn- We will reserve the symbol (g> for tensor products with 
respect to the A-B subdivision, and the symbol M for tensor products regarding 
the I-II subdivision. Strong superadditivity of the EoF is the inequality 

E F {p) >? E F (pi) + E F (pn), (7) 

for p a state on H, and pi and pn its reductions to Hi and Hn, respectively. 
The following Lemma is simple but crucial: 

Lemma 1. For any bounded function f defined on S(H), strong superadditivity 
off 

Kp)>U(pi) + f(pii), (8) 

is equivalent to subadditivity of the conjugate function f* with respect to the 
Kronecker sum: 

f*(x 1 mi + imx 2 ) <?f*(x 1 ) + r(x 2 ). (9) 

Proof. Set Z = X x IE I + I IE X 2 . Then, using © and assuming the validity of 
© yields 

f(p)= S u V Tr[ P X}-r(X) 
x 

> sup Tr[ P Z]-f*(Z) 

Xi .X2 

> sup Tr[piX 1 +puX 2 }-r(X 1 )-r(X 2 ) 
= Kpi) + Kpii), 

which is JHJ. 

The converse follows from ©. Assuming the validity of (JHJ yields 

f*(Z) = maxTr[pZ]-/(p) 
p 

< maxTr[yO/Xi + pnX 2 ] - f(pi) - f(pn) 
p 

= maxTr[piAi + p 2 X 2 \ - f( Pl ) - f(p 2 ) 

Pl,P2 

= nx) + f*(Y), 

which is □ 

The appearance of the Kronecker sum in Lemma 1 suggests that the consid- 
eration of the function /* o log is a more natural setting for studying additivity. 
Defining g := f* o log and setting Xj, = logMj, 1© becomes 

g{M x M M 2 ) <lg{Mi) + g(M 2 ), 
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for Mi,M 2 G B + (H). Restating © and © in terms of M, we have 

g(M)= max Tr[plog(M)]-/(p) (10) 

f(p) = max Ti[plog(M)]-g(M), (11) 

Strictly speaking, these quantities are defined only for positive M. However, 
when M is singular, we can still make sense out of it by the usual extension 
Tr[jolog(M)] = —oo for any p that is not completely supported on the range of 
M. 

We can now restate Lemma 1 in the form of a Theorem, which is our first 
main result: 

Theorem 1. For any junction f defined on S(TL), and with g defined on B + {TL) 
by M(]\) . strong superadditivity of the convex closure f, 

f(p)>U(px) + f(pn), (12) 

is equivalent to subadditivity of g, 

g{Mx m M 2 ) <? g{M x ) + g(M 2 ). (13) 

Note that the expression Tr[/?log(A/)] — g{M) is invariant under multiplication 
of M by a positive scalar. Hence, one could impose the restriction Tr M = 1, i.e. 
that M should be a state, or alternatively M < I, which is what we shall do. 

An immediate corollary of this Theorem is the equivalence of the strong 
superadditivity of the EoF with the subadditivity of g — E* o log, where E* is 
the conjugate of the entanglement functional E(p) = S^Tr^ p). We have chosen 
to present Theorem 1 in the more general way because it obviates the rather 
remarkable independence of the Theorem on any property of the function / at 
all. Specifically, while for the sake of defining the EoF it is necessary to split up 
the Hilbert space into two parties A and B, this is something the Theorem is 
completely oblivious of. 

The only interesting feature of E we can exploit at this level is its concavity. 
Concavity allows to simplify the conjugation expression by replacing the max- 
imisation over all mixed states by a maximisation over pure states. Indeed, the 
argument of the maximisation in 

g(M)= max Tr[/jlog(M)] - E(p) 
pes(H) 

is a convex function of p, and it is well-known 9 that a convex function achieves 
its maximum over a convex set always in an extreme point of that set, in this 
case in a pure state. Thus: 

g(M) = maxTr[^log(M)] - E{V). 

View 

Theorem 1 reduces the additivity problem for the convex closure, originally 
defined as a minimisation over ensembles, to an equivalent problem for the con- 
jugate function, defined as a maximisation over pure states. If counterexamples 
are found for (15} , this automatically disproves strong superadditivity Ijl2|l . so 
this simplification does not come at the cost of reduced power. Specifically, by 
"inverting" the proof of Lemma 1 (or Theorem 1) and employing Proposition 1, 
we easily get the following: 
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Proposition 2. If p violates strong superadditivity of f, M\ is optimal for 

pi in fiij) . and M% is optimal for pn, then M\ M M2 violates subadditivity of 
g HBjj. If Mi M M 2 violates fT^j and p is optimal for Mi H M 2 in then p 

violates fJjjj) . 



5. Maximal Output Purity 

Exploiting Theorem 1, we will now show that strong superadditivity of Ep would 
follow as a consequence of another additivity conjecture, concerning quantum 
channel capacities. Recollect that, since E is concave, the optimal p in (@J will 
be an extreme point of the feasible set, i.e. a pure state, so: 

E*(X) = maxTr[^X] - E(&). (14) 

From the additivity of E over pure states it easily follows that the correspond- 
ing function g — E* o log is superadditive, hence subadditivity of g implies its 
additivity. 



5.1. Step 1. The maximisation in g can be rewritten in terms of a maximal 
eigenvalue A max : 

Lemma 2. For any M £ B + (H), 

g{M) := maxCTr^logM] - S(Tr A &)) 

= max A max (logAf + log(I A ®T)). (15) 

t£S(H b ) 

Note that we will henceforth consider log M+log(lU®T) as an operator restricted 
to the range intersection ran(M) n ran(I eg) r). 
Proof. 

max A max (log M + log I4, ® r) 

r 

= maxmaxTr[^(logM + log Ia r)] (16) 

T 1p 

= max max Tr[^ log M] + Tr[Tr^(^) log r] 

r ip 

= maxTr^logM] - S(Tr A V). (17) 

i> 

In step l|16|l we have used the Rayleigh-Ritz representation of a maximal eigen- 
value, and in step l|17J) we have used the fact that relative entropy is non-negative 
and attains the value zero when (and only when) its arguments are equal. Specif- 
ically: 

= min 5(p||r) 

T 

= min —S(p) — Tr[plogr] 

r 

= —S(p) — max Tr [p log t] . 

T 

□ 
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5.2. Step 2. Using the Lie- Trotter formula, the logarithm can be replaced by a 
limit of a power function. 

Lemma 3. 

exp 5 (M) = lim&V*(M), 
p— >o e 

where 

h p {M) := max||AP /2 (I ® t) p M p/2 \\ 

T 

and \ \.\\ denotes the operator norm. 

Proof. Taking the exponential of both sides of (|15|l and noting expA max (M) = 
A max exp(M), we get 

exp g(M) = max 1 1 exp(log M + log(I ® r))| |. 

r 

To make sense of this formula, we extend exp(logM 4- log(I <g> t)) as on the 
complement of ran(M) n ran(I eg) t), as in J^]. The Lie- Trotter formula has a 
continuous version (see the remark after Lemma 3.3 in 

exp{A + B) = lim ( exp(pA/2) exp(pS) exp(pA/2)) 1/p . 

p— >o 

In particular, this gives us 

exp(logM + log(I <g> r)) = lim (M p/2 (I ® t) p M p/2 ) 1/p . (18) 

Define the shorthand functions 

/(r) := ||exp(logM + log(I®r))|| 
f p (T) := \ \{M p/2 (I®t) p M p/2 ) 1/p \\ 

over S(Tt). By (|18|) and the triangle inequality for norms, f p converges pointwise 
to /. The functions f p are clearly continuous for p > 0. By Lemma 4.1 of \Rj\, f 
is continuous too. From JSj (p. 118) we have that f p decreases monotonously to 
/ asp decreases to 0. The set S(H), over which / and f p are defined, is compact. 
Hence, all the prerequisites are fulfilled to apply Dini's theorem [TJ|, and we get 
that the convergence of f p to / is uniform over S(7i). 

Finally, uniform convergence is equivalent with convergence in the sup-norm. 
By the triangle inequality for norms, that in turn implies that the sup-norm of 
f p converges to the sup-norm of /. Therefore, hp^ p (M) = max T / p (r) = H/pHs 
converges to ||/||s = max T /(r) = expg(M). □ 

Additivity of g would thus follow as a consequence of multiplicativity of h p , 
h p {M 1 MM 2 ) =? h p (M l )h p (M 2 ), for p j 0. Following [TH], we say that a property 
holds for p I a if it holds for an arbitrarily small, but finite, interval p G (a, a + e], 
e>0. 
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5.3. Step 3. The quantity h p {M) is formally equal to the maximal output purity 
|18U19U2()| of quantum filtering operations. Indeed, 

h p (M) = maxTr[<P(M p/2 (I®T) p M p/2 )] 

T,(f> 

= maxTr[T p Tr A [M p/2 $M p/2 ]] 

T,(f> 

= max 1 1 Tr A [M pj2 <PM pj2 \ \ \ q 

where q = 1/(1 — p) and ||.|| g denotes the Schatten g-norm 2B> an d Vq(A) is 
the maximal output purity measured by the Schatten g-norm of the (non-trace 
preserving) operation 

A : p i-> A(p) = Tr A [M p ' 2 pM p / 2 ] . (19) 
If this operation would be trace preserving, we would call it a channel. 

5.4- Step 4- We now claim that there is no advantage in restricting attention 
to operations of the form (fH^I . It is of course true that the class of operations 
(|19fl is rather specific. They admit a Kraus representation such that the block 
column matrix (Ai)i obtained by stacking the Kraus element matrices Ai ver- 
tically, equals M p / 2 , which is a positive matrix. Necessary conditions are that 
A\Ai — M p (which is < I) and the input dimension of the operation should 
equal the output dimension times the number of elements. 

However, as regards the maximal output purity question, these structural 
peculiarities offer no additional mileage. To see this, consider the specific case 
that M is a partial isometry M = USW, where £ = |1)(1| ® Is and U is any 
unitary, then 

v q (A) = max\\Tr A [US p / 2 U^US p / 2 U^}\\ q 

= max ||Tr A [?7Z ,p / 2 <r£ p / 2 ?y t l|L 

pen 

= max ||Tx A [[/(|l)(l|®#")E/ + ]|| g) 

<p £n-B 

which is the generic case for operations from H. A to H. A . Thus, the case for 
the "special operations" TC >— > TL A contains the generic TC A i— > H A case and is 
therefore not easier to prove. 

5.5. Step 5. The exponent p of M, occurring in A, is coupled to q 7 occurring in 
v q , via the relation q = 1/(1 — p). To cap off our argument, we "decouple" p 
and q by replacing M p / 2 with a general matrix < X < I, strenghtening our 
multiplicativity conjecture ever so slightly. Noting finally that p J. corresponds 
to q I 1, we get our second main result: 

Theorem 2. If v q (A) is multiplicative for q J, 1 and for any filtering operation 
A, then the entanglement of formation is strongly subadditive. 
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Multiplicativity of v q had been conjectured in |18j for trace preserving chan- 
nels. It has been proven for entanglement breaking channels |2()j . unital qubit 
maps |22| and depolarising channels |23| . but, unfortunately, was refuted in |19j 
for q > 4.79. Nevertheless, the conjecture might still be true for q J, 1. 

Theorem 2 has to be compared to the main technical result in , which states 
that additivity of the Holevo capacity for given channels implies additivity of 
the EoF for certain states. In a sense, our Theorem 2 is stronger because we get 
the stronger outcome of strong subadditivity. On the other hand, this comes at 
the price of having to consider non-trace-preserving operations. 

After the appearance of the first draft of this manuscript, Shor proved [2"4"| 
the equivalence of four additivity conjectures: strong superadditivity of the EoF, 
ordinary additivity of the EoF, additivity of the maximal output purity vs of 
a channel as measured by the entropy and additivity of the classical (Holevo) 
capacity of a channel. As multiplicativity of v q (A) for q J, 1 implies additivity of 
v s(A) Shor's third equivalence provides an alternative proof for our result 
Theorem 2. 



6. Conclusion 

In conclusion, we have shown how a simple convex analytical argument leads 
to a simpler formulation of the entanglement of formation and an especially 
simple equivalent condition for strong superadditivity of the EoF. Based on this 
we have found the second result that strong superadditivity of the EoF would 
follow as a consequence of the multiplicativity of the maximum output purity 
v q of quantum filtering operations, for q J, 1. 

Acknowledgement. We gratefully acknowledge comments by M.B. Plenio, J. Eisert, M.B. 
Ruskai and Ch. King. SLB currently holds a Wolfson-Royal Society Research Merit Award. 
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